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ABSTRACT 

Aims. The goal is to test the structure of hot metallic stars, and in particular the structure of a near-surface convection zone using 
asteroseismic measurements. Indeed, stellar models including a detailed treatement of the radiative diffusion predict the existence of 
a near-surface convection zone in order to correctly reproduce the anomalies in surface abundances that are observed in Am stars. 
Methods. The Am star HD 209625 was observed with the Harps spectrograph mounted on the 3.6-m telescope at the ESO La Silla 
Observatory (Chile) during 9 nights in August 2005. This observing run allowed us to collect 1243 radial velocity (RV) measure- 
ments***, with a standard deviation of 1.35 ms"'. 

Results. The power spectrum associated with these RV measurements does not present any excess. Therefore, either the structure of 
the external layers of this star does not allow excitation of solar-like oscillations, or the amplitudes of the oscillations remain below 
20-30 cms"' (depending on their frequency range). 

Key words. Stars: individual: HD 209625 - Stars: interiors - Stars: oscillations - Binaries: spectroscopic - Techniques: radial 
velocities 



' 1. Introduction 

In the spectral range from A- to F-type stars, almost 70% of 
• the non-chemically peculiar (CP) stars are currently classified as 
' 6 Set variables, while most non-variable stars are Am-type stars. 
, Although observations do not completly rule out variability in 
Am stars (or A-type metallic-line stars), they do rigorously con- 
strain it: either variability is extremely rare, or the pulsations are 
of very low amplitude. 

At their arrival on the ZAMS, Am stars have a Hell convec- 
tion zone, typical of A-F stars. Due to their low rotational veloc- 
ity, mixing processes are not efficient enough to prevent He from 
settling down. Helium leaves the external layers and is therefore 
no longer able to excite pulsations typical for the classical in- 
stability strip. Moreover, selective radiative acceleration causes 
the accumulation of iron-peak elements where these elements 
dominate the opacity, and leads to the appearance of iron-peak 
convection zones centered at a temperature of approximately 
200000K (Turcotte et al.'2000^ When assuming that there is a 
sufficient overshoot to homogenize the surface regions between 
the H, He and iron-peak convection zones, the predicted surface 
abundance variations closely resemble the abundance anomalies 
observed in Am stars (Richer et al. I2000I I. Detailed models of 
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Am stars thus predict the existence of a surface convection zone 
which may extend deeply into the external layers of the star For 
instance, in a typical A-type star of L7M0, the surface convec- 
tion zone is predicted to only include the outer 10 to 10 Mq 
of the star, while it may include the outer 10"^ to 10"^ of the 
star for an Am star with iron-peak convection zones (Richard et 
al. 1200 11 1. The total mass of the convection zone for such an Am 
star is then found to be very similar to the mass of the surface 
convection zone of a I.4M0 star at solar metallicity which is 
able to drive solar-like oscillations as observed in the subgiant 
Procyon A (see e.g. Martic et al. 120041 and Eggenberger et al. 
I2004l l. 

In order to probe the structure of the predicted convection 
zone due to iron overabundances, we performed high-accuracy 
asteroseismic observations of an Am star. Indeed, depending on 
its structure, this convection zone may be able to drive solar-like 
oscillations similar to the ones detected for less massive stars 
(see e.g. the review by Bedding & Kieldsen l2006l l. 

Typical A-type stars show only a few spectral lines (as 
opposed to F or G stars) and rotate fast, which prevents us 
from obtaining high-precision radial velocity (RV) measure- 
ments. However, this is not the case for Am stars characterized 
by under-abundances of calcium and scandium, and by over- 
abundances of iron-group elements. The large number of metal- 
lic absorption lines in Am spectra is particularly valuable to ob- 
tain precise radial velocities, needed to detect solar-like oscilla- 
tions. Moreover, the rotational velocity distribution of Am stars 
is very diff'erent from that of regular A stars (Debernardi, 2002ll: 
their mean rotational velocity lies below 50kms"', thus some 
Am stars have very low rotational velocities, a necessary con- 
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dition to obtain accurate radial velocities. Therefore, Am stars 
are the only hot stars for which the high radial velocity accuracy 
needed to detect solar-like oscillations can be reached. 

The target selected is the Am SBl binary HD 209625 
(A5/A9/F2; with a V magnitude of 5.28). This object is repre- 
sentative of the Am main sequence stars while having a quite 
low rotational velocity vsin/ - 9.6 km s"' (deduced from our 
radial velocity measurements). 

In this paper, we report on Doppler observations of 
HD 209625 performed with the Harps spectrograph, and result- 
ing in the determination of a very high-accuracy time series for 
an Am star. 

The observations and the associated data reduction are pre- 
sented in Sect. |2l along with a study of the binarity of the target. 
The power spectrum analysis is discussed in Sect. [3] followed by 
the conclusions in Sect.|4] 

2. Observations and study of binarity 

The Am star HD 209625 was observed during 4 half-nights 
in September 2004 and during 9 nights in August 2005 with 
the Harps spectrograph (Pepe et al. 120021 see e.g. Carrier 
& Eggenberger '2006' for solar-like oscillations with this in- 
strument) mounted on the 3.6-m telescope at the La Silla 
Observatory (ESO, Chile). We took consecutive exposures, each 
lasting between 140 and 300 s (depending on the airmass and 
the extinction), and separated by dead times of 30 s. In total, 
430 and 1243 spectra were collected (for the first and second 
run respectively), with a typical signal-to-noise ratio (S/N) rang- 
ing from 100 to 300 at 530 nm. In parallel to the stellar ex- 
posures, the spectrum of a thorium lamp carried by a second 
fiber was simultaneously recorded, in order to monitor the spec- 
trograph's stability. The spectra obtained were extracted on- 
line. The radial velocities were computed online by weighted 
cross-correlation with a numerical mask constructed from a 
G2 dwarf spectrum. In addition, they were also determined by 
weighted cross-correlation with different masks built from more 
appropriate synthetic spectra (spectral type between A5 and F2, 
see Carrier et al. 2002) and by the optimum-weight procedure 
(Connes II985I Carrier et al. |2005al l, although without signifi- 
cant gain. The use of a G2 template can introduce a small sys- 
tematic shift in the radial velocities but without consequence for 
the variability analysis of this star. 

The Am star HD 209625 being a SBl binary, it is neces- 
sary to correct the observations for the Keplerian movement 
before studying the presence of p-mode oscillations. In order 
to accurately determine the period of this binary system, we 
used all Harps radial velocity measurements, along with pre- 
vious CoRAVEL measurements (Debernardi 120021 ) and with data 
from Abt and Levy ( 119851 ). Then, the other parameters were de- 
termined using only the Harps measurements, because they are 
far more accurate. The eccentricity was found to be very small 
(0.001) and was thus fixed to zero. The derived orbital parame- 
ters are listed in Table [T] The RV modulation due to the binarity 
of the system is shown in the bottom panel of Fig. [1] while the 
top panel shows the RV measurements during the second run (the 
9 nights in August 2005) with the orbital modulation removed. 

We note that most nights seem to exhibit the same shape: an 
increase of the RV, followed by a decrease. In fact, the radial ve- 
locities are correlated both with the airmass and with the ratio 
of the flux in the blue part of the spectrum to the flux in the red 
part. This modulation is related to the atmospheric extinction. 
At the accuracy of the ms"', all spectrum orders show system- 
atic differences in velocity. The time series computed from the 
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Fig. 1. Harps radial velocity measurements of HD 209625. The 
bottom panel shows the radial velocities due to the orbital mo- 
tion, with the parameters of Table [T] The top panel presents the 
RV measurements from which the sinusoidal fit was removed, 
and after the average of each night was set to zero. 

Table 1. Orbital elements derived for HD 209625 when taking 
into account Harps, Coravel and Abt and Levy measurements 
(see text for more details). 



Parameter 


Value 


P (days) 


7.83238 ± 0.00002 


To (HJD) 


53420.2304 ± 0.0001 


e 


(fixed) 


ID (deg) 


(fixed) 


Vo(km s-') 


20.2630 ± 0.3 


Ki (km s-i) 


7.2150 ±0.4 


fli sin i (Gm) 


0.77708 ± 0.00004 


/i(m)(10-3 Mo) 


0.30550 ± 0.00005 


o-(o-c) (m s-') 


2.1 



cross-correlations, which are weighted means by the flux of cor- 
relations of each order, are thus correlated to the flux change. 
No corrections were applied in the online pipeline to take into 
account such atmospheric extinction changes. However, this ef- 
fect should be corrected in the new online version of the pipeline 
(F. Pepe, private communication). 

Although part of the analysis is based on both observing 
runs, the second run measurements are mainly used, because 
the complicated window function of the first run consisting of 
4 half-nights leads to a less precise result, however in full agree- 
ment with the first run. 

A very small variation, related to either intrinsic variations of 
the star or perhaps to a tiny eccentricity, remains in the data (i.e. 
shift of 1-2 ms~' between nights). For this reason, the mean 
value for each night was also subtracted. This could slightly 
modify the power spectrum at very low frequencies, but has no 
influence in the range where oscillations are expected (Carrier et 
al. i2005bi ). 

3. Power spectrum analysis 

The power spectrum of the velocity time series was com- 
puted using the Lomb-Scargle modified algorithm (Lomb [19761 
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Fig. 2. Power spectrum of HD 209625 (right), with a zoom on the low-frequency part (left). The white line indicates the mean noise 
level, while dotted and dashed lines represent an amplitude of 2 cr and 3 cr respectively. 



Table 2. Distribution of the 1243 Doppler measurements during ^ 
the second run. i 



Date 


Nb of spectra 


Duration 
(hours) 


Dispersion a 

(m s-') 


2005/08/01 


134 


7.63 


2.09 


2005/08/02 


154 


8.01 


1.99 


2005/08/03 


169 


8.57 


1.94 


2005/08/04 


71 


7.46 


1.97 


2005/08/05 


72 


8.66 


2.27 


2005/08/06 


170 


8.77 


1.60 


2005/08/07 


167 


8.61 


2.67 


2005/08/08 


156 


7.90 


1.63 


2005/08/09 


150 


8.23 


3.05 
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Fig. 3. Power density of HD 209625 (with the power axis in log- 
arithmic scale). The dashed, dotted and full lines correspond re- 
spectively to the first, the second and both runs. 



Scargle 119821) . and is shown in Fig. |2l Its time scale gives a 
formal resolution of 1.39 //Hz. As typically observed in such a 
power spectrum, the noise is made of two distinct components: 

- At high frequencies the noise is flat, indicative of the Poisson 
statistics of photon noise. The mean white noise level crpo„ 
calculated between 1.5 and 2.5 mHz is 0.005 11 m^ s"^, cor- 
responding to (Tamp = V^^pow"*^i74 = 6.3 cms"' in am- 
pHtude (Kjeldsen & Bedding |T9931 l. With 1243 measure- 
ments, this high frequency noise corresponds to ctrv - 
^Ncr^^jA - 1.3 ms"'. This radial velocity uncertainty 
is exceptionally low for an A-star, but remains nevertheless 
over the photon noise uncertainty limit which is estimated by 
the pipeline at 0.5-1. ms ' . This small discrepancy could ei- 
ther be explained by an intrinsic variability of the star (such 
as granulation or stellar activity) or by an instrumental effect. 

- Towards the lowest frequencies, the power spectrum scales 
inversely with the square of the frequency, as expected for 
instrumental instabilities. 

The mean noise is indicated by the white line in Fig. |2] along 
with the 2 cr value (dotted line) and 3 cr value (dashed line). 

As discussed by Kjeldsen & Bedding ( 1995), the conversion 
to the power density requires multiplying the power spectrum by 
the effective observing time, calculated by integrating under the 



spectral window. Density spectra were built for both runs sep- 
arately and also for the combined data. The effective observing 
time has a value of 0.87, 2.50 and 3.56 d for the first run, the 
second one and the combined data, respectively. The smoothed 
power density spectra are displayed in Fig. [3] The density spec- 
trum of the first run is the most affected by noise because of the 
smaller number of measurements. 

Scaling from the solar case (Kjeldsen & Bedding 1995), the 
frequency of the greatest mode of HD 209625 would be expected 
to lie near 500-600 yuHz, although this value cannot be trusted 
as the internal structure of an Am star is totally diff'erent from 
the one of a solar-type star Solar-like oscillations are detected 
neither in the power spectrum (see Fig. ^ nor in a smoothed 
version of the density spectrum (see Fig.[3]l at frequencies below 
2500 /iHz. Moreover, Fig. |2]shows that only two peaks (near 400 
and 1100/iHz) have an amplitude greater than 3cr, a result in 
agreement with pure noise (Carrier et al. 2005a). 

Although no oscillations could be detected by searching for a 
power excess directly in the periodogram or in the power density 
spectrum, we also tried to identify a comb-like structure in the 
power spectrum corresponding to solar-like oscillations, and in 
particular to ^ = modes, by computing autocorrelations (region 
by region and with different thresholds). However, no autocorre- 
lation showed the presence of such a comb-like structure. 
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Fig. 4. Simulations giving the probabiltiy of detection of the os- 
cillations versus the mean amplitude of the three highest modes. 
The oscillations are detected if an excess appears above a certain 
threshold in the density spectrum. This threshold varies between 
3 cr (black circle) and 5 cr (black triangle), in steps of 0.5 cr (see 
text for more details on the simulations). 



In order to determine a tighter upper limit on the oscilla- 
tions, we also performed simulations: artificial modes were in- 
troduced while reproducing the overall noise level. The frequen- 
cies of these injected modes were separated by 20yuHz, while 
their amplitudes were described by a Gaussian enveloppe with 
a cr of 50//Hz; this corresponded to a few peaks over a small 
frequency range. The main free parameter in these simulations 
was the amplitude of the Gaussian. The other two parameters 
were also modified to check the independance of the results: no 
significant change was observed provided no more than 3 or 4 
introduced peaks had a non-negligible amplitude. A Gaussian 
noise, randomly drawn (using Monte-Carlo simulations) and ad- 
justed to give the same value as in the observed power spec- 
trum, was then added to the time series. Using the same time 
sampling as observed, radial velocities were determined by as- 
suming an infinite lifetime for the modes. Since the power spec- 
trum of solar-like oscillations is known to depend on the poorly- 
known lifetime of the modes, we thus worked with the density 
spectrum, which is far less sensitive to the lifetime value. The 
density spectrum is expected to be flat if there is no signal. The 
noise introduces some modulation, characterized by the cTrm^ of 
the density spectrum. For each value of the introduced Gaussian 
amplitude, 1 000 different simulations were performed, leading 
to 1 000 density spectra. Different thresholds were used to de- 
termine the detection probability, ranging from 3cri-ms to Scr^^^ 
by steps of O-Sctj^s. Then, the number of density spectra where 
peaks larger than a given threshold were observed provided the 
detection probability, for this given threshold and Gaussian am- 
plitude of the modes. Finally, the mean of the 3 highest ampli- 
tude modes was computed. Fig. |4] shows the detection probabil- 
ity of modes with a given mean amplitude, for different detection 
thresholds. In the observed density spectrum (Fig. O, no peak 
above 3.3cr„„5 is seen. This 3.3cr„„s threshold, when compared 



to the results of the simulations, means that oscillations with a 
mean amplitude of 3.05 (respectively 2.70) times the noise can 
be detected with a probability of 99 % (respectively 95 %). 

4. Conclusion 

We do not detect any solar-like oscillations in the Am star 
HD 209625: either this object does not exhibit any oscillations, 
or the amplitude of its oscillations is below our detection level. 

To unambiguously detect a power excess due to solar-like os- 
cillations (at the 99 % confidence level), the noise in the Fourier 
transform must be, at least, 3.05 times smaller than the amplitude 
of the oscillations. Therefore, the oscillation amplitude is be- 
low 20 cm s"' if the star oscillates at frequencies above 800 //Hz, 
and could reach 30 cms"' for oscillations with frequencies near 
200juHz. Such solar-like oscillations therefore remain, if they 
exist, very limited in amplitude. 

This non-detection of solar-like oscillations with a clear up- 
per limit puts strong constraints on the structure of near-surface 
convection zones predicted by theoretical models of Am stars. 
Indeed, these observations allow us to rule out all models of the 
Am star HD 209625 that predict a surface convection zone that 
would be able to excite solar-like oscillations with an amplitude 
larger than 30 cm s"' . 
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